SF  298  MASTER  COPY 


KEEP  THIS  COPY  FOR  REPRODUCTION  PURPOSES 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMBNO.  0704-0188 


Public  reporting  burden  for  this  coltectton  of  information  is  estimated  to  average  1  hour  per  response,  including  the  lime  for  reviewing  instructions,  searching  existing  data  sourcea, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comment  regarding  this  burden  estimates  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1215  Jefferson 
Davis  Highway.  Suite  1204,  Arlington.  vA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Papemvork  Reduction  Project  (0704-0188),  Washington.  DC  20503. 


5  the  Office  of  Management  and  Budget.  Papemvork  Reduction  Project  (0704-0188),  Washington.  DC  20503. 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

20  Feb  97 


3.  REPORT  TYPE  AND  DATES  COVERED 


4.  TITLE  AND  SUBTITLE 

The  Effect  of  Electromagnetic  Fields  on  the  Stability  of 
Cylindrical  Jets  and  Rods 


6.  AUTHOR(S) 

David  L.  Littlefield 


5  FUNDING  NUMBERS 


7.  PERFORMING  ORGANIZATION  NAMES(S)  AND  ADDRESS(ES) 

Southwest  Research  Institute 

PO  Drawer  28510 

San  Antonio,  TX  78228-0510 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-221 1 


11.  SUPPLEMENTARY  NOTES 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


ftlt6  .I0-E6- 


The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  be  construed  as 
an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 


12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 


12  b.  DISTRIBUTION  CODE 


Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Maximum  200  words) 


The  effect  of  electromagnetic  fields  on  the  stability  of  cylindrical  jets  and  rods  was  studied  in  this  research. 
The  goal  of  the  program  was  to  develop  a  better  understanding  of  the  disruptive  mechanisms  that 
electromagnetic  fields  might  have  on  shape-charge  jets  and  kinetic  energy  penetrators.  A  combined 
analytical  and  experimental  approach  was  undertaken.  Stability  analyses  were  first  performed  on  perfectly 
plastic,  uniformly  elongating  jets  when  subjected  to  axial  electric  currents,  and  later  extended  to  the  case 
of  elastic-plastic  cylindrical  rods.  Results  from  the  analyses  were  verified  using  a  series  of  small  scale, 
high  voltage  capacitive  discharge  experiments. 


iMTC  QtTALlTY  INSPECTED  2 


19970321  020 


14.  SUBJECT  TERMS  .  _  . 

Stability,  Magnetohydrodynamics,  Plasticity,  Linear  Perturbation 

Theory,  Shape-Charge  Jet,  Kinetic  Energy  Projectile 


15.  NUMBER  IF  PAGES 

19 


16.  PRICE  CODE 


17  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20,  LIMITATION  OF  ABSTRACT 
OR  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIHED  UNCLASSIFIED  UL 


NSN  7540-01-280-5500 


Standard  Form  298  (Rav.  2-89) 

Prescribed  by  ANSI  Std.  239-18 

298-102 


TABLE  OF  CONTENTS 


Page 

1.0  STATEMENT  OF  PROBLEM .  1 

2.0  SUMMARY  OF  IMPORTANT  RESULTS .  2 

2.1  Thermomechanical  and  Magnetohydrodynamic  Stability  of  Cylindrical  2 

Stretching  Jets . 

2.2  Magnetomechanical  Instabilities  in  Elastic-Plastic  Cylinders .  3 

2.3  Computational  and  Experimental  Studies  of  Magnetomechanical  Instabilities 

in  Cylindrical  Metal  Rods .  3 

2.3.1  Introduction .  3 

2.3.2  Mathematical  Model .  5 

2.3.3  Idealized  Motion .  6 

2.3.4  Perturbation  Equations .  6 

2.3.5  Solutions  to  Perturbation  Equations .  7 

2.3.6  Experiments . .;... .  7 

2.3.7  Comparisons  Between  Theory  and  Experiment .  12 

2.3.8  Conclusions .  13 

3.0  LIST  OF  PUBLICATIONS  AND  TECHNICAL  REPORTS  PUBLISHED .  14 

4.0  OTHER  REFERENCES .  15 

5.0  LISTING  OF  SCffiNTMC  PERSONNEL  SUPPORTED  ON  THE  PROJECT .  16 

David  L.  Littlefield .  16 

Donald  Grosch .  16 

Jerry  Nixon .  16 

RayBurgamy.. .  16 


11 


LIST  OF  FIGURES 

....  .  Page 

Figure  1 .  Normalized  perturbed  radius  vs.  normalized  time  f  at  various  wavenum¬ 
bers  for  App  =  100,  Q  =  0.88,  E  =  0.005 1  and  R  =  2. 1 .  2 

Figure  2.  Normalized  frequency  ©j  vs.  normalized  wavenumber  k  for  plastic  distur¬ 
bances  in  a  cylinder  with  O  =  1/2,  'P  =  0  and  A  =  4 .  4 

Figure  3.  Normalized  growth  rate  vs.  normalized  wavenumber  k  for  plastic  distur¬ 
bances  in  a  cylinder  with  4>  =  1/2  and  T^  =  0 .  4 

Figure  4.  Model  for  Stability  Calculations .  6 

Figure  5.  Normalized  Frequency  ©j  versus  Normalized  Wavenumber  k  for  Perturba¬ 
tions  in  an  Elastic-Plastic  Cylinder  with  O  =  1/2,  T  =  0.01  and  A^p  =  1 .  8 

Figure  6.  Normalized  Growth  Rate  ©r  versus  Normalized  Wavenumber  k  for  Pertur¬ 
bations  in  an  Elastic-Plastic  Cylinder  with  $  =  1/2  and  'P  =  0.01 .  8 

Figure  7.  Schematic  of  Experimental  Setup .  9 

Figure  8.  Voltage  and  Current  Profiles  from  Test  39 .  10 

Figure  9.  Digital  Image  (Displayed  in  Reverse  Video  for  Clarity)  Showing  the  Rod 

Disruption  in  Test  39,  58  iis  After  Current  Initiation .  10 

Figure  10.  Digital  image  of  the  rod  disruption  from  Test  47  at  13  p,s  (displayed  in 

reverse  video  for  clarity) .  11 

Figure  11.  Digital  image  of  the  rod  disruption  from  Test  48  at  15  jxs.  The  Image  is 

Partially  Obscured  by  the  Shadow  of  a  B-dot  Probe .  12 


Hi 


1.0  STATEMENT  OF  PROBLEM 


The  effect  of  electromagnetic  fields  on  the  stability  of  cylindrical  jets  and  rods  was  studied 
in  this  research.  The  goal  of  the  program  was  to  develop  a  better  understanding  of  the  disruptive 
mechanisms  that  electromagnetic  fields  might  have  on  shape-charge  jets  and  kinetic  energy 
penetrators.  A  combined  analytical  and  experimental  approach  was  undertaken.  Stability  analyses 
were  first  performed  on  perfectly  plastic,  uniformly  elongating  jets  when  subjected  to  axial  electric 
currents,  and  later  extended  to  the  case  of  elastic-plastic  cylindrical  rods.  Results  from  the  analyses 
were  verified  using  a  series  of  small  scale,  high  voltage  capacitive  discharge  experiments. 
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2.0  SUMMARY  OF  IMPORTANT  RESULTS 


2.1  Thermomechanical  and  Magnetohydrodynamic  Stability  of  Cylindrical  Stretching  Jets 

Stability  calculations  performed  previously  by  Littlefield  [4.1  -  4.6]  showed  that  electro¬ 
magnetic  fields  can  alter  the  stability  characteristics  of  isothermal,  uniformly  elongating  plastic 
jets.  The  analyses  demonstrated  that  axial  electric  currents  always  make  the  jet  more  unstable  and 
increase  the  growth  rate  of  naturally-occurring  instabilities,  even  when  small  currents  are  employed. 
These  results  were  then  expanded  to  include  high  levels  of  electric  current,  where  thermal  energy 
transfer  mechanisms  become  important  [3.1,  3.2].  It  was  discovered  that  thermal  softening  of  the 
jet  due  to  Joule  heating  has  atendency  to  increase  the  growth  rate  of  perturbations.  This  is  illustrated 
in  Figure  1,  where  a  plot  of  the  magnitude  of  perturbations  to  the  radius  of  the  jet  are  shown  as  a 
function  of  time  for  various  initial  wavelengths.  Perturbations  to  the  radius  are  normalized  by  their 
initial  value,  and  time  is  normalized  by  the  inverse  of  the  axial  strain  rate.  Values  of  dimensionless 
parameters  App,  Q,  E  and  R  were  100, 0.88, 0.005 1 ,  and  2. 1 ,  respectively.  Physically  App  represents 
a  ratio  of  electromagnetic  to  inertial  forces,  Q.  a  ratio  of  inertial  to  plastic  forces,  E  a  ratio  of  kinetic 
energy  to  sensible  heat  and  R  a  ratio  of  magnetic  convection  to  diffusion.  As  is  evident  from  the 
figure,  the  fastest  growing  disturbances  have  associated  wavenumbers  of  about  2n,  which  corre¬ 
sponds  to  a  wavelength  equal  to  one  half  the  jet  diameter.  A  detailed  description  of  the  analysis 
and  additional  results  can  be  found  in  refs.  [3.1,  3.2,  3.4]. 


Figure  1.  Normalized  perturbed  radius  f  vs.  normalized  time  f  at  various 
wavenumbers  for  App  =  100,  Q  =  0.88,  E  =  0.0051  and  R  =  2.1 
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2.2  Magnetomechanical  Instabilities  in  Elastic-Plastic  Cylinders 

Constitutive  models  used  in  the  stability  analyses  were  then  expanded  from  perfectly  plastic 
behavior  to  the  more  general  case  of  elastic-plastic  flow.  Results  from  the  models  indicated  that 
stability  of  elastic-plastic  rods  is  significantly  altered  when  axial  electric  currents  are  introduced. 
In  general,  the  frequency  of  oscillating  perturbations  is  suppressed,  and  the  growth  rate  of  distending 
instabilities  is  hastened.  Both  of  these  tendencies  appear  to  be  intensified  once  plastic  flow  has 
initiated.  This  is  illustrated  in  Figs.  2  and  3,  where  plots  of  the  normalized  frequency  (O;  and 
normalized  growth  rate  are  shown,  respectively,  versus  the  normalized  wavenumber  k.  The 
values  of  the  parameters  O  and  'F  were  0.5  and  0,  respectively.  In  Fig.  2,  the  parameter  A^p  was 
set  to  4,  and  in  Fig.  3,  growth  rates  for  a  range  of  values  for  A^p  is  shown.  Here  the  parameters  O, 
T,  and  Agp  physically  represent  a  ratio  of  elastic  to  compression  forces,  a  ratio  of  plastic  to  com¬ 
pression  forces,  and  a  ratio  of  electromagnetic  to  compression  forces,  respectively.  The  frequency 
spectrum  shown  in  Fig.  2  exhibits  substantial  suppression  when  compared  to  the  same  result  in  the 
absence  of  electric  current.  Furthermore,  the  growth  rates  in  Fig.  3  indicate  that  disturbances  which 
grow  in  amplitude  with  time  can  occur,  whereas  in  the  absence  of  current  distending  perturbations 
are  not  possible;  the  growth  rates  of  these  disturbances  and  the  wavenumber  at  which  the  maximum 
growth  rate  occurs  tends  to  increase  with  the  level  of  electric  current.  A  complete  description  of 
these  analyses  and  a  more  detailed  discussion  can  be  found  in  refs.  [3.3,  3.5,  3.6],  and  for  brevity 
will  not  be  repeated  here. 

23  Computational  and  Experimental  Studies  of  Magnetomechanical  Instabilities  in 

Cylindrical  Metal  Rods 

The  stability  analyses  of  elastic-plastic  rods  was  further  continued  to  consider  the  effects  of 
finite  electrical  conductivity,  where  it  was  also  found  that  introduction  of  the  current  can  cause  both 
fluctuating  and  distending  instabilities  in  the  rod.  A  comparison  of  the  results  to  similar  calculations 
for  perfectly  conducting  rods  predicted  that  the  effect  of  conductivity  is  to  increase  the  growth  rate 
is  distending  perturbations  and  lower  the  frequency  of  oscillatory  disturbances.  To  validate  the 
results  from  the  theory,  a  series  of  experiments  was  performed  by  discharging  electrical  energy 
from  a  high  voltage  capacitor  bank  through  thin  copper  rods.  The  development  of  axisymmetric 
distending  instabilities  in  the  rod  was  verified  using  a  series  of  high  speed,  digital  images.  Com¬ 
parisons  with  the  model  were  reasonable  and  within  the  limitations  of  the  analytical  theory.  These 
results  have  recently  been  reported  in  two  conferences  [3.7, 3.8],  but  have  not  yet  been  reported  in 
the  archival  literature,  and  for  convenience  are  presented  here  in  more  detail. 

23.1  Introduction 

The  instability  and  breakup  of  cylindrical  metal  rods  when  subjected  to  large  electromagnetic 
fields  has  become  and  area  of  increased  interest  over  the  past  several  years.  An  understanding  of 
this  phenomenon  is  critical  to  several  industrial  and  military  applications,  including  high  voltage 
circuit  interruption,  electromagnetic  forming,  magnetically-confined  fusion,  and  projectile  dis- 
raption.  For  example,  it  is  well  Imown  that  when  large  electric  currents  are  passed  through  cylindrical 
conductors  that  the  cylinder  has  a  tendency  to  explode,  breaking  up  into  a  large  number  of  thin, 
straited  segments  [4.7].  Although  rapid  expansion  of  the  plasma  formed  from  the  vaporized  wire 
certainly  contributes  to  the  observed  breakup,  electromagnetically-induced  instabilities  may  also 
be  partially  responsible  for  the  phenomenon.  The  same  processes  may  also  provide  a  means  for 
effective  protection  of  armored  vehicles;  for  example,  if  a  large  electric  current  was  introduced  into 
a  projectile,  it  might  be  disrupted  in  much  the  same  way  as  an  exploding  wire.  However,  in  order 
to  utilize  this  technology  effectively,  a  complete  understanding  of  the  coupling  of  electromagnetic 
effects  to  the  mechanical  behavior  of  the  material  must  be  obtained.  In  this  paper,  the  theory 
governing  the  formation  and  propagation  of  fundamental  magnetomechanical  instability  modes  in 
metal  cylinders  is  examined.  The  analysis  is  compared  to  results  from  small  scale  experiments 
performed  using  copper  rods. 
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Figure  3.  Normalized  growth  rate  (0^  vs.  normalized  wavenumber  k  for 
plastic  disturbances  in  a  cylinder  with  0  =  1/2  and  ^  =  0. 
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In  the  analysis,  the  stability  of  an  infinitely  long,  metal  cylinder  is  considered.  An  axial 
electric  current  is  introduced,  creating  an  azimuthal  magnetic  field  in  the  surrounding  vacuum,  and 
interacting  with  the  current  to  generate  an  electromagnetic  pinching  force,  directed  towards  the  axis 
of  the  cylinder.  The  conductivity  of  the  cylinder  is  assumed  to  be  small  enough  so  that  the  current 
remains  completely  diffused  throughout  the  cross  section;  this  condition  is  more  realistic  than  the 
assumption  of  perfect  conductivity  made  in  previous  analyses  [3.3, 3.5, 3.6].  The  electromagnetic 
force  gives  rise  to  a  condition  of  uniform  uniaxial  strain  in  the  cylinder.  Two  dimensional  axi- 
symmetric  perturbations  to  this  motion  are  considered,  and  the  equations  governing  the  motion  of 
these  disturbances  are  derived  using  linear  perturbation  theory.  These  equations  are  solved  to 
determine  how  the  amplitude  of  smdl  disturbances  changes  over  time. 

Results  from  a  series  of  small  scale  experiments  are  presented  to  validate  the  theoretical 
predictions.  Copper  wires  are  subjected  to  an  axial  electric  current  caused  by  discharge  from  a  bank 
of  high  voltage  capacitors.  Evolution  of  the  instability  and  breakup  of  the  cylinder  is  recorded 
through  a  series  of  high-speed,  digital  images. 


2.3.2  Mathematical  Model 

Consider  a  long,  cylindrical  metal  projectile  that  is  deformed  in  its  interior  as  a  result  of  an 
applied  electric  current.  The  geometry  of  the  cylinder  is  depicted  in  Fig.  4.  The  radius  at  an  arbitrary 
point  along  the  surface  of  the  rod  is  denoted  as  rj,(z,r),  and  the  unit  normal  vector  pointing  into  the 
surrounding  vacuum  is  n.  The  equations  governing  axisymmetric  motion  in  the  interior  of  the 
projectile  can  be  expressed  as 


9v,.  V,.  1  dp 

dz  K  dt’ 
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^ 

dr^  dr  dz^ 


(4) 


where  p  is  the  density;  K  is  the  bulk  modulus;  |i  is  the  majgnetic  permeability;  v,  and  are  the  radial 
and  axial  components  to  the  velocity,  respectively;  p  is  the  pressure;  Xy  is  a  component  of  the 
deviatoric  shear  stress  tensor;  and  Bq  is  the  azimuthal  component  of  the  magnetic  field.  Eq.  (1)  is 
the  conservation  of  mass,  Eqs.  (2)  and  (3)  the  conservation  of  radial  and  axial  momentum, 
respectively;  and  Eq.  (4)  the  magnetic  diffusion  equation.  When  coupled  with  the  appropriate 
constitutive  relations  in  the  cylinder,  Eqs.  ( 1 )  -  (4)  represent  a  complete  system  of  equations  to  solve 
for  the  motion  in  the  cylinder  and  in  the  surrounding  vacuum. 
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Figure  4.  Model  for  Stability  Calculations. 


2.3.3  Idealized  Motion 


The  governing  equations  given  in  the  previous  section  admit  a  solution  under  uniaxial  strain 
conditions.  This  idealized  motion  is  characterized  by  a  uniform  radial  compression  of  the  rod  as  a 
result  of  the  Lorentz  force  created  by  the  applied  electric  current.  Assuming  the  current  is  sufficient 
in  magnimde  to  initiate  yielding  through  the  entire  cross  section  of  the  cylinder,  the  pressure  and 
radial  displacement  become 


P(i  =  KpO--r^) 


'-f 


'PO 


(5) 


where  p  =  pIK,  =  yJ^IAi^rloK,  O  =  GIK,  ^  =  YIG,  G  is  the  shear  modulus,  Y  is  the 

yield  strength,  and  the  0  subscript  is  used  to  imply  a  zeroeth-order  solution.  Here  is  the  radial 
displacement  as  is  obtained  from  integrating  the  velocity  with  respect  to  time.  The  dimensionless 
parameters  that  appear  in  these  solutions  are  A^p,  O,  and  'P,  which  physically  represent  a  ratio  of 
electromagnetic  to  compression  forces,  elastic  shear  to  compression  forces,  and  plastic  to  elastic 
shear;  respectively. 


2.3.4  Perturbation  Equations 


Solutions  for  the  idealized  motion  in  the  rod  given  in  the  previous  section  are  now  subjected 
to  a  small  disturbance.  The  equations  governing  the  time  evolution  of  these  disturbances  are  derived 
by  using  linear  perturbation  Aeory.  Solutions  to  these  equations  are  used  to  identify  perturbations 
that  grow,  decay  or  oscillate  in  amplitude  with  time. 


The  perturbation  equations  are  formed  in  the  usual  way  by  expanding  the  independent  vari¬ 
ables  in  a  power  series  of  the  form/=/o  +/;£  +  ... ;  where  the  0,  1  and  2  subscripts  denote 

zeroeth,  first  and  second  order  approximations  to  the  variable/,  respectively;  then  retaining  only 
terms  which  are  of  order  £.  A  natural  consequence  of  this  expansion  is  that  the  resulting  linear 
differential  equations  can  be  Fourier  analyzed  with  respect  to  time  and  the  axial  dimension,  so  that 
each  of  the  Fourier  modes  depends  on  a  normalized  growth  rate  and  ®  and  wavenumber  k,  that  is 
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(6) 


/i=/(Oexp(cof+A:f), 


where  f  =  t/{rj,o'\JplK),  z  =  zirho,  i  =  V-  1  and /  is  the  amplitude  of /;.  After  the  expansions  are 
formed  and  the  equations  are  solved,  nontrivial  solutions  are  found  to  exist  only  when  a  certain 
dispersion  relation  is  satisfied;  this  relation  describes  the  (co,  k)  space  where  solutions  exist  to  this 
eigenvalue  problem.  The  details  of  the  derivation  and  the  resulting  dispersion  relation  are  described 
elsewhere  [3.1 1],  and  for  brevity  will  not  be  repeated  here. 

2.3.5  Solutions  to  Perturbation  Equations 

Solutions  to  the  dispersion  relations  described  in  the  previous  section  were  generated  for  <I> 
=  1/2  and  'F  =  0.01,  which  are  typical  values  for  many  metals.  The  dispersion  relations  also  admit 
solutions  for  real  and  imaginary  values  for  ©;  imaginary  values  correspond  to  oscillatory  instability 
modes  and  real  values  to  distending  modes.  In  the  discussion  that  follows,  ©j  is  the  imaginary 
component  to  ©  and  ©^  the  real  component,  and  are  referred  to  as  the  normalized  frequency  and  the 
normalized  growth  rate,  respectively. 

In  Fig.  5,  a  plot  of  the  stability  spectrum  for  oscillatory  modes  is  shown  when  A^p  =  1,  where 
©,  is  shown  as  a  function  of  the  wavenumber  k.  The  normalized  frequency  and  wavenumber  are 

related  to  the  dimensional  frequency /and  wavelength  as  f=  (A{sKlp  Httr^  and  =  Inr^lk, 

respectively.  Solutions  are  plotted  along  lines  that  exhibit  similar  deformation  characteristics,  and 
thus  belong  to  the  same  deformation  mode.  There  are  many  similarities  between  these  results  and 
the  results  shown  elsewhere  for  elastic  and  elastic-plastic  rods.  For  example,  the  frequency  of 
perturbations  tends  to  fluctuate  with  increasing  wavenumber,  as  is  seen  in  perfectly  conducting 
elastic-plastic  rods  [3.6],  but  there  is  an  overall  tendency  for  the  frequency  to  increase  with  k,  as  is 
observed  in  elastic  rods  [3.5,  3.9].  Suppression  of  the  frequency  occurs  for  the  lowest  oscillatory 
modes,  as  was  observed  for  the  perfectly  conducting  cylinder;  however,  the  magnitude  of  the 
suppression  is  significantly  larger.  For  example,  when  ^  =  10  the  normalized  frequency  of  the 
lowest  mode  shown  in  Fig.  5  is  1.105,  whereas  the  frequency  for  a  perfectly  conducting  rod  under 
the  same  conditions  is  3.961;  and  when  no  electric  current  is  present  this  frequency  becomes  4.041. 
The  higher  frequency  modes  are  also  influenced  by  the  conductivity;  albeit  to  a  much  lesser  degree 
than  the  lowest  mode.  For  example,  when  k  =  n/2  the  normalized  frequencies  for  the  second  mode 
are  2.000  and  2.132  in  the  present  case  and  for  perfect  conductivity,  respectively. 

The  introduction  of  electric  current  also  produces  distending  instabilities  in  the  rod.  Shown 
in  Fig.  6  is  the  normalized  growth  rate  ©,  versus  wavenumber  for  several  different  values  of  A^p. 

The  normalized  growth  rate  is  related  to  the  dimensional  growth  rate  a  as  a^dip^KIp/r^o.  It  is 
evident  from  these  results  that  the  perturbations  grow  in  amplitude  at  a  faster  rate  as  the  current  is 
increased.  The  fastest  growing  perturbation  has  a  wavelength  that  tends  to  decrease  with  increasing 
electric  current.  Similar  results  are  also  observed  in  perfectly  conducting  rods  [3.6];  however,  for 
any  given  Aep  the  magnitude  of  the  growth  rates  are  larger  by  a  factor  of  two  or  more  in  the  present 
case. 

2.3.6  Experiments 

In  order  to  investigate  the  implications  of  the  theory,  a  series  of  small  scale  experiments  was 
performed.  A  68  pF,  30  kV  capacitor  bank  was  used  to  discharge  large  electric  currents  through 
small  copper  rods.  A  schematic  of  the  major  components  in  the  experimental  setup  are  shown  in 
Fig.  7.  The  capacitor  bank  is  connected  across  a  set  of  aluminum  Walker  plates  [4.8]  using  a  set 
of  four  RG-220  coaxial  cables,  and  a  1.02  mm  diameter,  OFHC  copper  test  cylinder  was  placed 
across  these  plates.  In  order  to  simplify  the  diagnostics,  the  tests  were  performed  quasistatically 
and  the  rod  remained  stationary  between  the  plates  for  the  duration  of  the  experiment.  However, 
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Figure  5.  Normalized  Frequency  cOj  versus  Normalized  Wavenumber  k  for 
Perturbations  in  an  Elastic-Plastic  Cylinder  with  <5  =  1/2,  T  =  0.01  and  A^p  =  1. 


Figure  6.  Normalized  Growth  Rate  co^  versus  Normalized  Wavenumber  k 
for  Perturbations  in  an  Elastic-Plastic  Cylinder  with  =  1/2  and  'P  =  0.01. 
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Figure  7.  Schematic  of  Experimental  Setup. 

this  does  not  eliminate  the  applicability  of  the  results  to  moving  projectiles,  since  the  electromagnetic 
forces  remain  unchanged  in  a  frame  of  reference  moving  along  with  the  rod.  The  electric  current 
was  initiated  using  a  25  kV,  100  kA  ignitron  connected  between  the  plates  and  ground. 

Several  forms  of  diagnostic  measurements  were  used  to  interpret  results  from  the  experiments. 
High  voltage  dividers  recorded  the  voltage  across  the  capacitor  bank  and  the  Walker  plates,  and 
Pearson  coils  to  measure  the  current  flowing  through  the  rod.  B-dot  probes  were  positioned  radially 
around  the  circumference  of  the  rod  to  measure  the  magnetic  field.  An  Imacon  468  high  speed, 
digital  camera  with  a  six-shot  capability  was  used  to  image  the  experiment.  The  rod  was  illuminated 
by  diffusing  the  beams  from  two  20  mW,  He-Ne  lasers,  and  a  632-nm  bandpass  filter  was  placed 
over  the  camera  lens  to  filter  out  all  wavelengths  other  than  the  light  produced  by  the  laser.  This 
scheme  helped  eliminate  much  of  the  "flash"  that  typically  occurs  in  high-voltage  discharge 
experiments. 

Voltage  and  current  profiles  produced  from  a  typical  low  voltage  experiment  are  shown  in 
Fig.  8,  where  the  current  and  voltage  across  the  Walker  plates  is  shown  versus  time  for  Test  39. 
The  initial  voltage  across  the  plates  is  about  7.2  kV.  Current  and  voltage  profiles  are  representative 
of  what  is  typicdly  seen  in  the  discharge  of  an  RLC  circuit.  The  peak  current  experienced  by  the 
rod  is  42.4  kA,  and  occurs  17  fis  after  initiation.  From  the  shape  of  the  decay  in  the  profile  beyond 
17  |is,  the  inductance  and  resistance  of  the  test  setup  can  be  estimated.  The  inductance  is  about 
850  nH,  but  the  resistance  tends  to  increase  with  time;  at  early  times  the  resistance  is  about  28  mQ, 
but  later  in  the  experiment  the  resistance  increases  to  roughly  45  mO.  A  similar  behavior  was  also 
observed  when  the  experiment  was  repeated  at  higher  voltages,  but  the  time  where  the  transition 
occurs  in  the  resistance  becomes  progressively  shorter. 

Disruption  of  the  rod  is  consistent  with  what  is  observed  in  the  voltage  and  current  profiles. 
Shown  in  Fig.  9  is  an  image  from  Test  39,  showing  the  perturbation  and  breakup  of  the  cylinder  58 
p,s  after  initiation  of  the  current.  Fig.  9  is  a  typical  image  showing  the  early  stages  of  disruption 
for  the  cylinder  at  this  particular  voltage.  Not  only  is  the  rod  is  beginning  to  break  up  into  small 
segments,  but  it  also  shows  indications  of  a  kink  instability  developing,  which  tends  to  bend  the  rod 
out  of  plane.  The  wavelength  of  the  segments  is  varied  and  falls  within  a  range  of  roughly  1  to  2 
times  the  radius  of  the  rod,  corresponding  to  wavenumbers  of  %  and  27c,respectively.  Arcing  has 
occurred  where  a  section  of  the  rod  has  vaporized.  Six  microseconds  later,  the  cylinder  almost 
completely  vaporizes,  with  indications  of  a  vaporous  halo  surrounding  remnants  of  the  original  rod. 
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Voltage  and  current  from  test  39 


Figure  9.  Digital  Image  (Displayed  in  Reverse  Video  for  Clarity)  Showing 
the  Rod  Disruption  in  Test  39, 58  ps  After  Current  Initiation. 
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The  change  in  resistance  inferred  from  the  current  profiles  also  occurs  at  the  same  time  the  rod 
appears  to  be  vaporizing;  thus,  it  seems  reasonable  to  assume  that  the  increase  is  caused  by  the 
copper  rod  transitioning  from  a  solid  and/or  liquid  state  and  into  a  plasma. 

When  the  current  is  increased,  characteristics  of  the  breakup  of  the  cylinder  change  dra¬ 
matically,  and  begin  to  look  much  like  what  is  typically  observed  in  exploding  wire  experiments 
[4.7].  Shown  in  Figs.  10  and  11  is  a  sequence  of  images  from  two  different  tests  depicting  the 
dismption  in  the  rod,  when  a  voltage  of  13.9  kV  is  applied  across  the  plates.  Fig.  10  is  from  Test 
47  at  13  )a,s;  and  Fig.  1 1  from  Test  48  at  1 5  ps.  During  die  early  stages  of  dismption,  the  main  effect 
of  the  current  is  to  increase  the  diameter  of  the  cylinder;  after  1 1  ps,  the  diameter  increased  to  3.49 
mm.  By  13  ps,  axisymmetric  perturbations  are  beginning  to  develop  in  the  interior  of  the  rod,  as 
shown  in  Fig.  1 0.  The  wavenumber  associated  with  these  perturbations  is  0.427 based  on  the  original 
radius  of  the  cylinder,  and  1 .47  based  on  the  radius  at  1 1  ps.  Fig.  1 1  indicates  that  these  perturbations 
grow  very  rapidly  in  amplitude  over  the  next  two  microseconds,  which  eventually  disperse  into  a 
plasma  cloud  by  17  ps. 


Figure  10.  Digital  image  of  the  rod  disruption  from  Test  47  at  13  ps 
(displayed  in  reverse  video  for  clarity). 
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Figure  11.  Digital  image  of  the  rod  disruption  from  Test  48  at  15  |xs.  The  Image 
is  Partially  Obscured  by  the  Shadow  of  a  B-dot  Probe. 

2.3.7  Comparisons  Between  Theory  and  Experiment 

Results  from  the  experiments  can  be  used  to  make  comparisons  with  the  predictions  from  the 
perturbation  theory.  The  analytical  model  assumes  the  rod  is  isothermal,  so  some  adjustments  to 
account  for  the  temperature  changes  need  to  be  applied.  A  simple  calculation  can  be  performed  by 
integrating  an  energy  equation  of  the  form 

— 

Bt  pcp\%rl 

where  T  is  the  temperature,  X  the  resistivity,  and  the  specific  heat.  Using  temperature-dependent 
values  for  X  and  c^,  together  with  the  current  profiles  from  the  experiments,  the  temperature  can  be 
determined.  Assuming  an  initial  temperature  of  300K,  the  temperatures  of  the  rod  become  505K 
at  1 7  jis  for  the  1 3.9  kV  experiment,  and  7 1  OK  at  60  |i.s  for  the  7.2  kV  test,  which  are  the  approximate 
times  in  each  experiment  when  disruption  in  the  rod  becomes  the  most  prevalent.  The  temperature 
increase  results  in  thermal  softening  of  the  rod,  which  can  be  accounted  for  in  the  model  by  making 
the  appropriate  adjustments  in  the  yield  strength.  Using  the  corrected  values  together  with  the 
current  profile  from  the  experiments,  eigenvalues  that  satisfy  the  dispersion  relations  governing  the 
stability  of  the  rod  [3.10]  can  be  determined.  Applying  the  conditions  from  the  higher  voltage  tests, 
these  dispersion  relations  predict  that  maximum  in  occurs  for  (©,,  k)  =  (0.1 186,  1.277).  This 


wavenumber  falls  within  the  range  of  values  measured  from  Tests  47  -  48.  However,  using  the 
conditions  from  Test  39  the  predicted  maximum  in  occurs  for  k)  =  (0.0346,  0.8595).  This 
wavenunfibef  is  not  in  agreement  with  the  values  seen  in  the  experiment,  where  values  of  3  to  6 
were  typically  observed. 

The  predicted  values  for  O),  can  also  be  compared  with  the  growth  rates  in  perturbations  seen 
in  the  experiments.  For  example,  the  amplitudes  of  the  disturbances  shown  in  Figs.  7a  -  7b  can  be 
measured  directly  and  correspond  to  a  growth  rate  of  approximately  5.4  x  10^  sec  \  On  the  other 
hand,  the  growth  rate  corresponding  to  co,  =  0.1 186  as  predicted  from  the  model  is  8.7  x  10^  sec  ', 
which  is  reasonably  close  to  the  measured  value.  The  growth  rate  from  Test  39  could  not  be 
measured,  but  seemed  to  be  quite  small  up  until  the  point  where  the  rod  appeared  to  be  vaporizing. 
The  predicted  value  using  the  model  is  2.5  x  10^  sec  ',  which  is  large  enough  to  be  measured  in  a 
test;  thus  the  model  overpredicts  the  growth  rate. 

It  seems  likely  that  limitations  of  the  theoretical  model  are  responsible  for  disagreements 
between  the  experimental  and  analytical  results.  For  example,  in  the  model  the  rod  is  assumed  to 
be  solid  and  isothermal.  Although  the  temperature  changes  predicted  from  the  experiments  seem 
low  enough  to  validate  this  assumption,  localization  of  the  current  from  asymmetries  or  skin  effects 
could  cause  "hot  spots",  resulting  in  localized  melting  and/or  vaporization.  Certainly  the  dynamics 
of  the  perturbation  growth  will  change  if  this  occurs.  Furthermore,  since  the  temperature  prediction 
at  breakup  was  higher  for  the  7.2  kV  experiments,  it  seems  likely  that  those  tests  might  be  more 
susceptible  to  localized  phase  changes;  this  may  explain  why  inadequate  agreement  was  obtained. 

2.3.8  Conclusions 

The  stability  of  elastic-plastic  cylinders  has  been  considered  in  this  analysis.  In  general,  for 
a  given  level  of  electric  current,  isothermal  conducting  rods  exhibit  both  distending  and  oscillatory 
perturbations.  The  growth  rates  of  distending  perturbations  tend  to  be  larger,  and  the  frequency  of 
oscillatory  disturbances  smaller,  than  is  obtained  for  perfectly  conducting  rods.  The  formation  of 
circulatory  motion  is  also  more  aggressively  suppressed  when  finite  conductivity  is  taken  into 
account. 

Comparison  of  the  analytical  predictions  with  experimental  results  is  reasonable,  suggesting 
that  instabilities  developing  while  the  rod  is  still  a  solid  may  be  responsible  for  the  observed  breakup. 
The  largest  disagreements  are  obtained  when  temperature  variations  in  the  rod  are  expected  to  be 
the  most  important.  Thus,  a  more  detailed  consideration  of  thermal  effects  is  needed  to  improve 
the  predictive  capability  of  the  model. 
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